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We compute exactly the asymptotic distribution of scaled height in a (l+l)-dimensional 
anisotropic ballistic deposition model by mapping it to the Ulam problem of finding the longest 
nondecreasing subsequence in a random sequence of integers. Using the known results for the Ulam 
problem, we show that the scaled height in our model has the Tracy- Widom distribution appearing 
in the theory of random matrices near the edges of the spectrum. Our result supports the hypothesis 
that various growth models in (1 + 1) dimensions that belong to the Kardar-Parisi- Zhang univer- 
sality class perhaps all share the same universal Tracy- Widom distribution for the suitably scaled 
height variables. 
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Growth processes are ubiquitous in nature. The past 
few decades have seen an extensive research on a wide 
variety of both discrete and contiuous growth models 
HI ■ A large class of these growth models such as the 



Eden model restricted solid on solid (RSOS) mod- 
els 0, directed polymers polynuclear growth models 
(PNG) [|| and ballistic deposition models (BD) 0] are be- 
lieved to belong to the same universality class as that of 
the Kardar-Parisi-Zhang (KPZ) equation describing the 
growth of interface fluctuations ||. This universality is, 
however, somewhat restricted in the sense that it refers 
only to the width or the second moment of the height 
fluctuations characterized by two independent exponents 
(the growth exponent /3 and the dynamical exponent z) 
and the associated scaling function. Moreover, even this 
restricted universality is established mostly numerically. 
Only in very few special discrete (1+1 )-D models, the 
exponents /3 = 1/3 and z = 3/2 can be computed ex- 
actly via the Bethe ansatz technique A natural and 
important question is whether this universality can be ex- 
tended beyond the second moment of height fluctuations. 
For example, does the full distribution of the height fluc- 
tuations (suitably scaled) is universal, i.e. is the same 
for different growth models belonging to the KPZ class? 
Moreover, the KPZ-type equations are usually attributed 
to models with small gradients in the height profile and 
the question whether the models with large gradients be- 
long to the KPZ universality class is still open. 

Recently Prahofer and Spohn found an exact 
mapping between a specific PNG model and the so- 
called longest increasing subsequence (LIS) problem, also 
known as the Ulam problem. The LIS problem was first 
raised by Ulam in the early 60's [ll|, then the interest 
in it reappeared in the mathematical literature in 70 's 
since the work of Vershik and Kerov [l2j. The exact 
mapping of PNG to LIS and the subsequent utilization 



of the exact results available for the LIS problem al- 
lowed Prahofer and Spohn to find (besides the exact KPZ 
growth exponent (3 — 1/3) the exact asymptotic height 
distribution in the PNG model [ljj. This distribution 
turned to be the well known Tracy- Widom distribution 
appea ring in the theory of edge states of random ma- 
trices |l3j . Around the same time, Johansson showed 
rigorously 0] that a specific (1+1)-D directed polymer 
model, believed to be in the KPZ universality class, also 
has the same Tracy- Widom distribution for the scaled 
height (energy) variable. Gravner et. al. found the same 
Tracy- Widom distribution in another class of (1 + 1) di- 
mensional growth models which they called 'oriented dig- 
ital boiling' model|l5|. It would be interesting to know 
whether there are other growth models such as the RSOS 
or the BD ones, which are believed to be in the KPZ uni- 
versality class as far as the second moment is concerned, 
also share the same Tracy- Widom distribution for the 
scaled height. 

The purpose of this Letter is to present a BD model 
which can be mapped exactly to the LIS problem and 
hence it shares the same Tracy- Widom distribution as 
the PNG model. This exact result, in combination with 
the results of [ToL [Til Hj| . then lends support to the hy- 
pothesis that perhaps all these different growth models, 
at least in (1+1) dimensions, share the same universal 
Tracy- Widom distribution for scaled height. This hy- 
pothesis, if true, puts the universality on a much stronger 
footing going beyond just the second moment. Inciden- 
tally, to our knowledge, our model is the first exact so- 
lution for the full asymptotic height distribution of BD 
type systems. 

Before describing our model, it is worth summarizing 
the main results for the LIS problem that we use later. 
Take a set of n integers {1, 2,3,..., n}. Consider all the 
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n! possible permutations of this sequence. For any given 
permutation, let us find out all possible increasing subse- 
quences (terms of a subsequence need not necessarily be 
consecutive elements) and from them find out the longest 
one. For example, take n = 10 and consider a par- 
ticular permutation {8,2,7,1,3,4,10,6,9,5}. From this 
sequence, one can form several increasing subsequences 
such as {8, 10}, {2, 3, 4, 10}, {1, 3, 4, 10} etc. The longest 
one of all such subsequences is either {1,3,4,6,9} as 
shown by the underscores or {2,3,4,6,9}. The length 
l n of the LIS (in our example l n = 5) is a random vari- 
able as it varies from one permutation to another. In 
the Ulam problem one considers all the nl permutations 
to be equally likely. Given this uniform measure over 
the space of permutations, what is the statistics of the 
random variable l n ? Ulam found numerically for the av- 
erage length (l n ) the asymptotic behavior (l n ) ~ c^/n for 
large n. Later this result was established rigorously by 
Hammersley |l6j and the constant c = 2 was found by 
Vershik and Kerov Recently, in a seminal paper, 

Baik, Deift and Johansson (BDJ) derived the full 
distribution of l n for large n. In particular, they showed 
that l n — > 2y/n + n x / & x f° r large n, where the random 
variable x has an n-independent distribution which hap- 
pens to be the Tracy- Widom distribution for the largest 
eigenvalue of a random matrix drawn from the Gaussian 
Unitary Ensemble ^3| . They also showed that when the 
sequence length n itself is a random variable drawn from 
a Poisson distribution with mean (n) = A, the length of 
the LIS converges for large A to 

? A ->2VA + A 1 / 6 x, (1) 

where x has the Tracy- Widom distribution. The fixed 
n and the fixed A ensembles are like the canonical and 
the grand canonical ensembles in statistical mechanics. 
The detailed form of the Tracy- Widom distribution is 
rather complicated and not very illuminating (see [To| 
for a picture). The BDJ results led to an avalanche of 
subsequent mathematical works . 

In our (1+1)-D BD model columnar growth occurs se- 
quentially on a linear substrate consisting of L columns 
with free boundary conditions. The time t is discrete and 
is increased by 1 with every deposition event. One starts 
at t = with an empty substrate. At any stage of the 
growth, a column (say the fc-th column) is chosen at ran- 
dom with probability p — and a "brick" is deposited 
there which increases the height of this column by one 
unit, H k — > Hf. + 1. Once this "brick" is deposited, it 
screens all the sites at the same level in all the columns to 
its right from future deposition, i.e. the heights at all the 
columns to the right of the fc-th column must be strictly 
greater than or equal to H k + 1 at all subsequent times. 
For example, in Fig^ the first brick (denoted by 1) gets 
deposited at t = 1 in the 4-th column and it immediately 
screens all the sites to its right. Then the second brick 
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FIG. 1: Growth of a heap with asymmetric long-range inter- 
action. The numbers inside cells show the times at which the 
blocks are added to the heap. 

(denoted by 2) gets deposited at t = 2 again in the same 
4-th column whose height now becomes 2 and thus the 
heights of all the columns to the right of the 4-th column 
must be > 2 at all subsequent times and so on. Formally 
such growth is implemented by the update rule, 

H k (t + l)=ma X {H k (t),H k _ 1 (t),...,H 1 (t)} + l. (2) 

The model is anisotropic and evidently even the average 
height profile (H k {t)) depends nontrivially on both the 
column number k and time t. Our goal is to compute the 
asymptotic height distribution P k (H, t) for large t. 

It is easy to find the height distribution Pi(H,t) of 
the first column, since the height there does not depend 
on any other column. At any stage, the height in the 
first column either increases by one unit with probability 
p = \ (if this column is selected for deposit) or stays the 
same with probability 1 — p. Thus P\ (H, t) is simply the 
binomial distribution, Pi(H,t) = (^)p h (l ~ pY~ H with 
H < t. The average height of the first column thus in- 
creases as (Hi(t)) — pt for all t and its variance is given 
by g\ (£) = tp(l — p). While the first column is thus 
trivial, the dynamics of heights in other columns is non- 
trivial due to the right-handed infinite range interactions 
between the columns. For convenience, we subsequently 
measure the height of any other column with respect to 
the first one. Namely by height h k {t) we mean the height 
difference between the (fc+l)-th column and the first one, 
h k (t) = H k+1 (t) - Hi(t), so that h (t) = for all t. 

To make progress for columns k > 0, we first consider 
a (2+l)-D construction of the heap as shown in FigEl 
by adding an extra dimension indicating the time t. In 
Fig|21 the x axis denotes the column number, the y axis 
stands for the time t and the z axis is the height h. In 
this figure, every time a new block is added, it "wets" 
all the sites at the same level to its "east" (along the 
x axis) and to its "north" (along the time axis). Here 
"wetting" means "screening" from further deposition at 
those sites at the same level. This (2 + 1)-D system of 
"terraces" is in one-to-one correspondence with the (1 + 
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column number, k 



FIG. 2: (2 + 1) dimensional "terraces" corresponding to the 
growth of a heap in FigQ 

1)-D heap in Fig^ This construction is reminiscent of 
the 3D anisotropic directed percolation (ADP) problem 
studied by Rajesh and Dhar [l^. Note however, that 
unlike the ADP problem, in our case each row labelled 
by t can contain only one deposition event [2^. 

The next step is to consider the projection onto the 
2D (x, y)-plane of the level lines separating the adjacent 
terraces whose heights differ by 1. In this projection, 
some of the level lines may overlap partially on the plane. 
To avoid the overlap for better visual purposes, we make 
a shift (x, y) — ► (x+h(x, y), y) and represent these shifted 
directed lines on the 2D plane in FigEl 

The black dots in Fig|31 denote the points where the 
deposition events took place and the integer next to a 
dot denotes the time of this event. Note that each row 
in Fig|3 contains a single black dot, i.e. only one de- 
position per unit of time can occur. In FigEl there are 
8 such events whose deposition times form the sequence 
{1,2,3,4,5,6,7,8} of length n = 8. Now let us read 
the deposition times of the dots sequentially, but now 
column by column and vertically from top to bottom 
in each column, starting from the leftmost one. Then 
this sequence reads {8,3,5,1,2,6,4,7} which is just a 
permutation of the original sequence {1, 2, 3, 4, 5, 6, 7, 8}. 
In the permuted sequence {8,3,5,1,2,6,4,7} there are 
3 LIS's: {3, 5, 6, 7}, {1, 2, 6, 7} and {1, 2, 4, 7}, all of the 
same length l n = 4. There is a greedy algorithm called 
the "patience sorting" game devised by Aldous and Dia- 
conis to determine this length of the LIS 01- This game 
goes as follows: start forming piles with the numbers in 
the permuted sequence starting with the first element 
which is 8 in our example. So, the number 8 forms the 
base of the first pile. The next element, if less than 8, 
goes on top of 8. If not, it forms the base of a new pile. 
One follows a greedy algorithm: for any new element of 
the sequence, check all the top numbers on the existing 
piles starting from the first pile and if the new number 
is less than the top number of an already existing pile, 
it goes on top of that pile. If the new number is larger 
than all the top numbers of the existing piles, this new 
number forms the base of a new pile. Thus in our exam- 
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FIG. 3: The directed lines are the level lines separating adja- 
cent terraces with height diffrence 1 in Fig. 2, projected onto 
the (x,y) plane and shifted by (x,y) — > (x + h{x,y),y) to 
avoid partial overlap. The black dots denote the deposition 
events. The numbers next to the dots denote the times of 
those deposition events. 

pie, we form 4 distinct piles: [{8, 3, 1}, {5, 2}, {6, 4}, {7}]. 
The number of piles (4) is same as the length l n = 4 of 
the LIS of this permuted sequence. In fact, Aldous and 
Diaconis proved that the length of the LIS l n is ex- 
actly equal to the number of piles in the corresponding 
patience sorting game. 

Let us note one immediate fact from FigEl The 
numbers belonging to the different level lines in 
FigCH are in one-to-one correspondence with the piles 
[{8, 3, 1}, {5, 2}, {6, 4}, {7}] in Aldous-Diaconis patience 
sorting game. Hence, each pile can be identified with 
an unique level line. Now, the height h(x,t) at any 
given point (x, t) in FigO is equal to the number of 
level lines inside the rectangle bounded by the corners: 
[0,0], [x, 0], [0,i], [x,t]. Thus, we have the correspon- 
dence: height = number of level lines = number of piles 
= length /„ of the LIS. However, to compute we need 
to know n which is the number of black dots inside this 
rectangle. 

Once the problem is reduced to finding the number 
of black dots or deposition events, we no longer need 
the Fig|3| (as it may confuse due to the visual shift 
(x, y) — > (x + h(x, y), y)) and can go back to Fig0 where 
the north-to-east corners play the same role as the black 
dots in Fig[3 In Fig[3 to determine the height h^(t) of 
the /c-th column at time t, we need to know the number 
of deposition events inside the 2D plane rectangle Rt.t 
bounded by the four corners [0,0], [k,Q], [0, t], [k,t]. Let 
us begin with the last column k — L. For k — L the num- 
ber of deposition events n in the rectangle Rh,t is equal 
to the time t because there is only one deposition event 
per time. In our example n = t = 8. For a general k < L 
the number of deposition events n inside the rectangle 
Rk t t is a random variable, since some of the rows inside 
the rectangle may not contain a north-to-east corner or 
a deposition event. The probability distribution Pk,t(n) 
(for a given [k,t]) of this random variable can, however, 
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be easily found as follows. At each step of deposition, a 
column is chosen at random from any of the L columns. 
Thus, the probability that a north-to-east corner will fall 
on the segment of line [0, k) (where k < L) is equal to 
k/L. The deposition events are completely independent 
of each other, indicating the absence of correlations be- 
tween different rows labelled by t in Fig0 So, we are 
asking the question: given t rows, what is the probabil- 
ity that n of them will contain a north-to-east corner? 
This is simply given by the binomial distribution 



Pk,t(n) 



(3) 



where n < t. Now we are reduced to the following prob- 
lem: given a sequence of integers of length n (where n 
itself is random and is taken from the distribution in 
Eq.®), what is the length of the LIS? Recall that this 
length is precisely the height hk(t) of the fc-th column at 
time t in our model. In the thermodynamic limit L — > oo 
for t ^> 1 and any fixed k such that the quotient A = ^ 
remains fixed but is arbitrary, the distribution in Eq.Q 



becomes a Poisson distribution P(n) 



-XX 



^r, with the 



mean A = ^. We can then directly use the BDJ result 
in Eq.Q to predict our main result for the height in the 
BD model, 



h k (t) 




(4) 



for large A = tk/L, where the random variable \ nas the 
Tracy- Widom distribution. Using the known exact value 
(x) = —1.7711... from the Tracy- Widom distribution^^ > 
we find exactly the asymptotic average height profile in 
the BD model, 



(**(*)> 



1.7721... — 



1/6 



(5) 



The leading square root dependence of the profile on the 
column number k has been seen numerically j^. The 
Eq.jSJ) also predicts an exact sub- leading term with A; 1 / 6 
dependence. Similarly, for the variance, cr£(t) = ([hk(t) — 

1/3 

(h k (t))} 2 ), we find asymptotically: vl(t) — > cq (^) , 
where c = ([x - (x)] 2 ) = 0.8132... [13]. Eliminating 
the t dependence for large t between the average and 
the variance, we get, cr\{t) ~ a(hk(t)) 213 where the con- 
stant a = co/2 2 / 3 = 0.51228... and = 1/3, thus re- 
covering the KPZ scaling exponent. In addition to the 
BD model with infinite range right-handed interaction 
reported here, we have also analyzed the model (analyt- 
ically within a mean field theory and numerically) when 
the right-handed interaction is short ranged j^. Surpris- 
ingly, we found that the asymptotic avera ge h eight profile 
is independent of the range of interaction 21] . 

In summary, we have shown that the asymptotic scaled 
height in an anisotropic (1+1)D BD model has the Tracy- 
Widom distribution. Our exact result, in combination 



with those of Refs. [ToL IlH Il5| where the same distribu- 
tion was found in rather different growth models, sug- 
gests that the universality in all these growth processes 
is perhaps much wider than it was thought before, ex- 
tending to the full asymptotic height distribution and is 
not restricted only to the second moment of the height. 

We thank the hospitality of the Institut Henri Poincare 
in Paris where this work was initiated during the 
trimester " Geometry and Statistics of Random Growth" 
in 2003. 
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